SOME MEASURE THEORY ON STACKS OF 
NETWORKS 

JACK MORAVA 

Abstract. By counting symmetries of graphs carefully (or equiv- 
alently, by regarding moduli spaces of graphs as zero-dimensional 
orbifolds), certain measures on these collections (elsewhere called 
'exponential random graphs') can be reinterpreted, with the aid of 
special cases of Wick's theorem, as Feynman-style measures on the 
real line. Analytic properties of the latter measures can then be 
studied in terms of phase transitions - in particular, in models for 
spaces of scale-free trees. 



Introduction: This paper is an attempt at an essentially elementary 
account of a single example, which provides evidence for the emergence 
of an interesting statistical mechanics of networks. 

For more than fifty years physicists have elaborated techniques in- 
troduced by Feynman, which organize calculations of probabilities on 
spaces of particle histories as sums over graphs. More recently [2], re- 
searchers have begun to use these methods in reverse, to attack prob- 
lems in enumerative combinatorics by analytical means. 

Here [§2.3] such (relatively familiar) ideas are formulated as an analog 
of the Fourier transform, as a way of evaluating sums over certain spaces 
of graphs (known in the literature [15] as exponential random measures) 
as integrals over the real line, with respect to certain Feynman/Gibbs- 
style measures. Collections of graphs or networks seem intuitively very 
discrete, but the transformed sums have interesting (and accessible) 
analytic properties, which can sometimes be described in ways similar 
to the phase transitions studied in thermodynamics. 

The first section below reviews standard background material, but it 
contains a significant technical point: sums over graphs are really in- 
tegrals over zero- dimensional moduli spaces, and if we treat orbifold 
points respectfully (ie by weighting their symmetries appropriately, by 
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what physicists call Faddeev - Popov determinants) we get better for- 
mulas. Similar issues appear in other apparently discrete contexts; for 
example, in Siegel's mass formula for quadratic lattices [18 VII §6.5 
(remark)]. 

Without compelling examples this would be empty formalism. Section 
three sketches the basic properties of a model (developed by physicists 
[3,4,12] interested in (among other things) the statistics of long-chain 
polymers) for a phase change in an ensemble of scale-free trees (ie 
having nodes of valence n occurring with power-law probability). This 
model has some remarkable features: its phase transition is of high (ie 
fourth) order, and its critical exponents depend to some extent on the 
parameters of the model. I hope this note will help make it, and its 
many generalizations, accessible for the closer study it deserves. 

I'd like to acknowledge helpful correspondence and discussions with 
Zdzislaw Burda, Pawel Gajer, Soren Galatius, Florin Spinu, and Ben 
Mann. Special thanks are due to Charles Epstein, who very kindly 
straightened out my considerable confusion about the polylogarithm 
[6]. 

This work was supported by DARPA's program on the fundamental 
questions of biology. 



I. The stack of weighted graphs 
This section assembles standard facts and notation from graph theory: 

1.1.0 The very familiarity of graphs can lead to confusion. I will be 
concerned with finite, abstract graphs; such a thing can be defined [8] 
to be a finite set G with an involution a, together with a retraction t : 
G — > G a onto the fixed point set of the involution. Elements of G a are 
the vertices of the graph, while the elements of the complement G — G a 
are said to be its half- edges; the quotient of this set by the involution 
is the set of edges. The valence of a vertex v G G a is the cardinality 
k(v) of t _1 (f ), ie the number of its incident edges; I will always assume 
this is positive. A vertex of valence one defines a terminal edge, or leg, 
of the graph; the remaining elements of G a comprise its set Ver°(G) of 
internal nodes. (Legs will always be terminal.) 
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These definitions are somewhat counter-intuitive but they have some 
virtue; a good test case is the problem of enumerating the possible 
labellings of a graph with one edge and two vertices. 

Graphs define a category, the morphisms being maps of sets with invo- 
lution, compatible with the specified retractions; but isomorphisms will 
play a particularly important role in this note. Abstract graphs have 
natural geometrical realizations as one- dimensional simplicial complexes; 
we will not use them here, but they are important in generalizations in 
which the graphs map to interesting configuration spaces [5]. 

1.1.1 I will always assume that an abstract graph G carries a specified 
order, or labeling, on its set of legs. I will also consider graphs weighted 
by a function wq from the internal nodes of G to {0, 1, 2, . . . }. I will 
write Aut(G) for the group of isomorphisms of G, compatible with its 
labeling and weight function (which will in general be suppressed from 
the notation). T, (resp. WT,) will denote the stack (or groupoid, or 
zero-dimensional orbifold) of abstract (resp. weighted) graphs, subject 
to a stability condition on internal nodes (made precise below), with 
isomorphisms as maps. |r.| and |Wr.| will be the associated sets of 
isomorphism classes of objects. 

1.1.2 It is often easiest for technical purposes to work with connected 
graphs, and I will write T [resp WT] for the subcategories of such 
things. The genus, or first Betti number, of a connected graph is g{G) = 
E — V + 1, with E its total number of edges, and V the number of 
vertices; this is additive under disjoint union, but for some purposes 
the Euler characteristic x{G) (= 1 — g{G) for connected graphs) which 
is also additive, is more useful. 

The number v{G) of (external) legs, and the number e = E — v of inter- 
nal edges, are other useful additive functions. The Euler characteristic 
(or genus), together with the number of legs, define a kind of bigrading 
on the category of graphs. More generally, the generalization 

g w (G):= J2 ™g(v)+ 9 (G) 

ueVcr°(G) 

of the genus is additive on weighted graphs, and defines a similar bi- 
grading on WT,. 

A connected G is stable with respect to its weight function, if 2(wg(v) — 
1) + k(v) > for each internal node [10 §2]; note that for unweighted 
graphs (ie with wq = 0) this precludes nodes of valence two. I'll write 
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|r(x, v) | for the set of isomorphism classes of connected graphs of Euler 
characteristic x with v legs, and so forth. 

There are many other interesting naturally defined functions on these 
sets, such as the degree sequence 

D(G) = {v = di,d 2 , ■ ■ ■ ,d m } 

which assigns to G, the partition of 2E with d k equal to the number of 
vertices of valence k. The distribution of weights, or numbers of legs, 
define similar functions on |WT(x, v)\. 

The function which sends G, weighted or not, to its group Aut(G) of 
automorphisms is another important example, as is the related group 
of automorphisms which are allowed to change the labelings on the 
legs. When G is connected, this is just the product of Aut(G) with 
the symmetric group of permutations of the labels, but if not, things 
are more complicated; see below. Note that (labeled) trees have no 
automorphisms. 

1.2 Graphs in general can be described in terms of the symmetric 
product of collections of connected graphs, but we will need to keep 
careful track of symmetries. 

1.2.0 An ordered n-tuple of elements from a set X is an element of the 
n-fold Cartesian product X n of X with itself. Permuting the order of 
these elements defines an action 

S"xl^ X n 

of the symmetric group on this set, and the quotient X n /E n is the set 
of unordered n-tuples of elements from X. The coproduct, or disjoint 
union, 

]JX n /E n := SP°°(X) 

n>0 

of these quotients is the free abelian monoid generated by X: its ele- 
ments can be written as formal finite sums 

in which the n» are non-negative integers, and the Xi are (not necessarily 
distinct) elements of X; of course these formal sums are subject to 
various rules, such as n{x}+{x} = (n+l){x}, etc. SP°°(X) is naturally 
graded by the degree ^2rii{xi} ^J2 n i- 

1.2.1 The isomorphism class of a general element G of T. can thus be 
written as a formal sum ^njGj], defined by the disjoint union of 
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copies of (the isomorphism classes of) connected graphs G{. In other 
words, 

|r.| = sp°°|r| 

as graded sets, with the understanding that we now admit a 'vacuum' 
graph with no vertices, cf eg [4 Ch 7]. If, in the presentation of G as 
a sum, the indexed components are distinct, then its automorphism 
group 

Aut(G) H^njAntiGi)) . 

is a kind of wreath product. If v = Y n i v i = Yl ^ r k, where r k is the 
number of components of G with precisely k terminal legs, then the 
group of automorphisms of G which are allowed to permute those legs 
will be a semidirect product of the restricted automorphism group with 

UK k - ni write 



m(G, v) = 



for the number of ways of labeling the legs of G. 

1.2.2 Finally, it may be worth noting (since the symmetric product 
construction is not much used in analysis) that a measure on X pushes 
forward to define a measure on SP°°(X). A measure on connected 
(weighted) graphs thus extends naturally to define a measure on general 
(weighted) graphs. 

The 'measures' of most interest in this paper will, however, take values 
in some field R((/c)) of formal series; k will play the role of a param- 
eter like Planck's constant in some asymptotic expansion. In fact the 
formulas will usually include even more parameters. 



2. A KIND OF FEYNMAN TRANSFORM 



This note is concerned with measures on spaces of (weighted) graphs, 
which is a subject with a large and somewhat disconnected literature; 
in [15 §5, 16] the term 'exponential random graphs' is used for a class 
of examples somewhat wider than those which will be in focus here. 

2.1 That general class of models employs a family ej, of functions on 
(isomorphism classes of) graphs, together with corresponding parame- 
ters Pk, called 'inverse temperatures'; a graph G is then assigned prob- 
ability 

P0(G) = exp(-£fte fc (G)) 
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(perhaps suitably normalized). If the e's are additive (under disjoint 
union of graphs), these probabilities will be multiplicative, and their 
values on general graphs will be extended from their values on con- 
nected graphs as described in the preceding section. 

That will be the case here, with some modifications: Vp will be defined 
by functions e w ^{G) which count the total number of vertices in G with 
weight w and valence k > 1, and I will also introduce parameters t, A 
and k, to keep track of the number of legs, edges and Euler character- 
istic (which are all of course additive). The literature of exponential 
random graphs allows more complicated e's, which count more general 
subconfigurations (various kinds of polygons, etc) of G; such data can 
be incorporated into measures of the form F{G)T > p{G). 

In this paper the key difference involves the enumeration of graphs with 
symmetries. If we define 

MG) = ^^exp(-^/3^(G)) 

then our 'measure' /i^a (actually taking values in the formal power 
series ring R((«))[[£, A]]) will be defined by summing the function 

G~Vp(G)^\^^ 

over subsets of |WT.|. The extra combinatorial factors in the definition 
of V/3 divide by the order of the group of extended automorphisms of 
the graph (allowed to permute the labels on the legs). Note that the 
order |Aut(G)| is not quite multiplicative under disjoint union. 

2.2 The interest of this class of measures is that under certain circum- 
stances they can be calculated, or at least described, in terms of formal 
measures defined by local interactions in a one-dimensional Euclidean 
space. 

Feynman introduced certain expressions of the general form 



/ 



interpreted as perturbed Gaussian integrals over infinite-dimensional 
spaces, and he showed that these could be in some sense evaluated 
in terms of formal series involving more and more complicated but 
finite-dimensional integrals; but it was the physicist Gian-Carlo Wick 
who systematized the combinatorics behind those calculations. This 
section is concerned with a drastically simplified case of Wick's result, 
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which provides a rigorous asymptotic expansion for expressions like the 
integral above. In this account the domain of integration will be taken 
to be one-dimensional [9, 10], but the techniques below generalize quite 
naturally to 'non-linear sigma models', in which the graphs are mapped 
by the parameter x to some kind of configuration space [7]. 

2.2.1 In our situation 

C(x) = \\- l x 2 - V(x) 

will be an analog of the Lagrangian function of physics. Its first term 
represents the quadratic term in the exponent of the background Gauss- 
ian integral, while 

fc>0 

plays the role of an 'interaction potential'; but we will allow coefficients 

which are power series in the asymptotic parameter. However, we need 
to require that 6 ,o — &o,i = ^0,2 = &i,o — 0, ie that the coefficients b g ^ = 
unless 2(g - 1) + n > (or, alternately: 1^(0) = mod k 2 , V'(0) = 
V"(0) = mod k); this signals a secret connection with moduli spaces 
of Riemann surfaces [9, 13]. This extra generality will not be used in 
the main example in §3, but it illustrates some of the flexibility of this 
class of models. 

2.2.2 Let Z(0) denote J R exp(— k^C^x)) dx (or, more precisely, its 
representation as an asymptotic series); then we can define the 'expec- 
tation value' of a function / on the line to be 

(/) ~ Z(O)- 1 / f{x) expi-K^Cix)) dx . 

JR 

If we now write 

f3 g ,k = - log b 9)k , 
then Wick's theorem takes the form of an asymptotic expansion 

(exp^te)) ~ e2^[l + J2 V ^' X{G)X£iG) ^G)^ 

Ge\wr+\ 

(independent of the domain of integration, as long as it contains the 
origin). The sum is taken over the set |W^r^"| of weighted 'non- vacuum' 
graphs: connected or not, but such that each component has at least 
one leg. From now on the leading term in right-hand expression will 
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be included in this sum, on the grounds that it can be regarded as the 
contribution of the empty set, understood as a non- vacuum graph. 

Note that when V — 0, this formula reduces to the classical fact that 
the Fourier-Laplace transform of a Gaussian function is again Gaussian. 

2.2.3 We can use some formal version 

(/>= / /(0<exp(«x)> 
Jr 

of Plancherel's theorem to extend this result. Since the asymptotic 
expression depends only on the germ of / at the origin, we will regard 
/ G as a power series; substituting t = in^ gives 

\WT+\ 

If 7fc(x) = x k /k\ is the fcth divided power of x, for example, then the 
coefficient 

<9^exp(|/tA<9 2 ) n/k(x)\ x =o = 7m(|«A) 
if k — v — 2m is even, and is zero otherwise. Note that the resulting 
"measures' are really formal sums of delta-functions and their deriva- 
tives, all with support at the origin! 

2.2.4 As formulated in [10], Wick's theorem in one dimension asserts 
that 

(27tk\)- 1/2 J expK- 1 (tx-\\- 1 x 2 +V(x^j dx ~ exp /sT 1 (|At 2 +V(t)^ 
where the formal transform 

has coefficients 

b XtU = |Aut(G)| _1 A e(G) J] b w(v)Mv) . 

Ge|r( x ,^)| u6Ver°(G) 

The term |At 2 will sometimes be interpreted as the contribution to V 
from the unique graph in T(l,2) (whose coefficient is not defined by 
the prescription above). 

The expansion in §2.2.2 just rewrites this: when t — 0, 

(2itk\)- 1/2 [ exp(-/t- 1 £(a;)) dx ~ exp(K- 1 1/(0)) , 
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with V(0) = Yl g >i bg,o K9 (since 2(g — 1) + n > 0). This is a sum over 
'vacuum graphs' (without legs), so the 'renormalized' sum V+(t) = 
V(t) — V(0) contains contributions only from graphs with at least one 
leg. Thus 

(exp(K-Hx)) = Z(0)~ 1 J expK- l {tx-C{x)) dx ~ exp K~\\\t 2 +V+(t)) , 

and the formula in §2.2 comes from expanding the exponential of 
k~ 1 V + : we get 

TT „* , in V- TT [Pa(G i )K-^X< Gi H v ^] ki 

n exp(iA Ut |-P, K -A-) = En 1 ;;.,,!,^ 1 

|r(x,^)|,^>i 

which, everything being more or less additive, rearranges into 

(where G — J] fcjGj, with the distinct). 

2.3 We can summarize this elementary calculation as follows: 
We consider two formal measure spaces, 

• the real line R, with respect to the formal measure 

He ~ (27mA)~ 1/2 / exp(-/t' 1 £(a;)) dx , 
Jr 

(supported at 0, as noted above), and 

• |Wr.|, with respect to the measure f/,p t \ defined by summing 

Let X v be the characteristic function of the set of non-vacuum weighted 
graphs with exactly v legs, and define 

m>0 

With this notation, we have the 
Theorem: The linear operator 

maps the space L 1 (R, fic) of formal functions isometrically to the formal 
span of the X v in L 1 (|VFr.|, pi/3,\). 
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This correspondence sees functions on the line only in the neighbor- 
hood of zero, and their transforms lie in the subspace spanned by the 
characteristic functions X v ; but we can nevertheless probe the space of 
graphs in more detail, by varying the parameters A and (5. 



3. Ensembles of scale-free trees 

This general framework seems very rich, even after rather drastic spe- 
cialization. This section is an account of some measures on spaces of 
trees defined by potentials 

^) = E^TT (= IWz) - * - 2-*- l x 2 ) 

fc>3 

of polylogarithmic form. [The 'missing' linear and quadratic terms are 
tracked by the parameters t and A -1 ; the conventions used here seem 
to simplify book-keeping in the long run [9].] The parameter n does 
not appear in this formula, so we will be concerned with unweighted 
graphs; in fact we will be mostly concerned with trees (ie graphs of 
genus zero, which however will not be rooted). 

This subject has a rather extensive literature in condensed-matter 
physics, closely related to the study of 'scale-free' networks [3,4,12], 
but one of the purposes of this note is to argue that it deserves more 
attention from mathematicians. In particular, it seems to me that this 
model is potentially at least as rich as (for example) the Ising model. 
It may also be useful to note that related path integral techniques have 
found other interesting applications in mathematical biology, eg [17]. 

3.1 The partition function for weighted graphs is the power-series ex- 
pansion of an analytic function 

Z(t,X) = {2tvk\)- 1/2 I exp(-K -1 (fa; - C(x))) dx 

near t = A = 0. Following ideas pioneered by Landau [1, 14 Ch VI], 
we can try to interpret its singularities in terms of phase transitions. 

The leading term in the statistical free energy for this ensemble of 
networks is the tree approximation 

.t" 



W(t) = limKlogZ ~ ±At 2 + lim V{t) = V Pp{G)\ e - 

Ge|r(i,i/)|,i/>2 



SOME MEASURE THEORY ON STACKS OF NETWORKS 



11 



the last equality uses the fact that labeled trees have no automor- 
phisms. We can omit the factorial in the denominator of the last term 
by writing it as a sum over isomorphism classes of unlabeled trees. 

3.2 A formal analog of Laplace's method for evaluating integrals with 
large parameter (ancestral to the saddle point or 'stationary phase' 
approximation), at a critical point x of 

<f>(x) = <f>(x ) + \<P"(xq) ■ (x - X ) 2 + . . . , 

is an expression of the form 

(2™)- 1 / 2 / exp(-K-V(aO) dx = 6XP( " K " 1 ^° )) [1+higher order in re ... ] . 
This implies a formula 

W{t, A) = tx - C(x ) 

for the tree level free energy, interpreted as a function of t, A via the 
equation t = C'(x ) for a critical point of tx — C(x). Rewriting that 
equation in iterable form yields 

x (t, A) = A(t + V'(x )) = At + XV'(Xt) + • • • G R[[t, A]] ; 

more generally, the implicit function theorem guarantees the existence 
of a solution x = xo(t,X) off the critical locus A -1 = V"(x ), ie away 
from 

t = x V"(x ) - V'(x ) ■ 

3.3 However, the polylogarithm has a subtle singularity at x — 1: 

Li Q+ i(e- x ) = T(-a)(x~iQ) a 

modulo smooth functions, where 

( x ~ i0) a ■= 1 [( x + i0) a + ( x - i0) a ] =x% + cos na x a _ 

is the average of two possible Gel'fand-Shilov regularizations of the 
complex power. [The first draft of this paper contained a very confusing 
sign error here; I'm very sorry about that. A fuller discussion of this 
formula has been since posted at [6 Prop. 3].] Thus, if the potential 
has the form specified at the beginning of this section, then at x — 1 
we recover equation (26) 

to = V"{1) - V\l) = 1 + C(a - 1) - 2C(a) 
of [3] for the locus of singularities in the model. 

As a varies, this defines an approximate hyperbola in the positive quad- 
rant of the (t, a) plane, partitioning it into two regions, as in Fig. 6 
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of [4]. Left of this line, the resulting measure is expected to be con- 
centrated on trees with an exponential degree distribution, while to 
its right it is thought to be supported on 'clumpy' concentrations, with 
many edges abutting on a relatively small number of hubs. The critical 
line defines a phase transition, with the measure concentrated on the 
scale-free trees. 

The behavior of function of A near the critical line varies, 

depending on whether a is greater or less than two. If a > 3, the first 
coefficient 

™=- x -i(t + V'(x ))+x^V"(x ) 

OXq 

in the Taylor series for A = X(x ) goes to zero near the critical line, so 

x = l-C(\ -\) 1 ' 2 + --- . 
On the other hand the saddle-point equation 

X- 1 = x 1 (t + V'(x )) 

can be rewritten 

A -1 = Ao 1 + x^(t + V'(x )) - (t + V\l)) , 
(with Ao = V'(l)); as t — > to this specializes to 

A- 1 = A 1 + (x^ - 1)V'(1) + x?(V'(x ) - V'{1)) . 
When 3 > a > 2, 

A" 1 ~ A 1 + T(l - a)(l - xqY' 1 + • • • , 

so 

x ~ 1 - C(A - A) 1/(q_1) H . 

It is tempting to take 

^logxo-tAo-A) 1 ^- 1 ) 
(when a G (2, 3)) as an order parameter for this transition [1,19 I §3]. 

Some further results about the higher-genus case can be found in [12]; 
I hope to return to these, and other questions about the generalized 
thermodynamical properties of these models, soon. 
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